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ABSTRACT 

Combining a new, increased dataset of 8 QSOs covering the Lya forest at redshifts 
1.5 < z < 3.6 from VLT/UVES observations with previously published results, we 
have investigated the properties of the Lya forest at 1.5 < z < 4. With the 6 QSOs 
covering the Lya forest at 1.5 < z < 2.5, we have extended previous studies in this 
redshift range. In particular, we have concentrated on the evolution of the line number 
density and the clustering of the Lya forest at z < 2.5, where the Lya forest starts to 
' show some inhomogeneity from sightline to sightline. We have fitted Voigt profiles to 

CO . the Lya absorption lines as in previous studies, and have, for two QSOs with z cm ~ 2.4, 

CN| ■ fitted Lya and higher order of Lyman lines down to 3050 A simultaneously. This latter 

approach has been taken in order to study the Ly/3 forest at z ~ 2.2 and the higher 
H 1 column density Lya forest in the Ly/3 forest region. 

For a given TVhi range, the Lya forest at 1.5 < z < 4 shows the monotonic 
evolution, which is governed mainly by the Hubble expansion at this redshift range. 
^ ■ In general, the Lya forest line number density (dn/dz) is best approximated with 

dn/dz = 6.1(1 + z )2- 4 7±o.i8 f or the Hi column density N B i = 10 13 - 64 ~ 17 cm~ 2 at 
1.5 < z < 4. When the results at < z < 1.5 from HS T observations are combined, 
5.^ ■ the slow-down in the number density evolution occurs at z < 1.5. For higher column 

^ I density clouds at 2Vki > 10 14 cm~ 2 , there is a variation in the line number density from 

03 sightline to sightline at z < 2.5. This variation is stronger for higher column density 

systems, probably due to more gravitationally evolved structures at lower z. The mean 
Hi opacity r H i is t H i(» = 0.0032 (1 + z )3-37±o.20 at 1.5 < z < 4. jjST observations 
^ . show evidence for slower evolution of thi at z < 1. For iV Hl = 1Q 12 - 5 ~ 15 cm~ 2 , the 

%— i ' differential column density distribution function, /(iVm), can be best fit by /(-/Vhi) oc 

Njii with ^ Ri 1.5 for 1.5 < z < 4. When combined with HST observations, the 
exponent (3 increases as z decreases at < z < 4 for iVm = 10 13 ~ 17 cm -2 . The 
correlation strength of the step optical depth correlation function shows the strong 
evolution from < z > = 3.3 to < z > = 2.1, although there is a large scatter along 
different sightlincs. The analyses of the Ly/3 forest at z ~ 2.2 are, in general, in good 
agreement with those of the Lya forest. 
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1 INTRODUCTION 



* The data used in this study are based on public data re- 
leased from the UVES Commissioning and Science Verification 
and from the OPC program 65.O-0296A (P.I. S. D'Odorico) at 
the VLT/Kueyen telescope, ESO, Paranal, Chile, 
t The line lists in the study are only available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5). 



The redshift evolution of the Lya forest imprinted in the 
spectra of high-z QSOs provides a powerful tool to probe 
the distribution and evolution of baryonic matter, and hence 
the formation and evolution of galaxies and the large scale 
structure, over a wide range of redshifts up to z ~ 6 (Sargent 
et al. 1980; Schaye et al. 1999; Kim, Cristiani & D'Odorico 
2001, hereafter KCD). 
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The evolution of the Lya forest is mainly governed by 
two physical processes. One is the Hubble expansion and 
another is the ionizing ultraviolet background flux (Theuns, 
Leonard & Efstathiou 1998; Dave et al. 1999; Schaye et al. 
2000; Bianchi, Cristiani & Kim 2001). At higher z, the Hub- 
ble expansion and the non-decreasing ultraviolet background 
cause a rapid evolution of the line number density per unit 
redshift, dn/dz. At lower z, the number of photons available 
to ionize the Lya forest becomes smaller, due to the decrease 
of the number of QSOs at z < 2, which is generally assumed 
to be the main source of the ionizing photons. As a result 
the rate of change of dn/dz with redshift is smaller. 

Studies on the forest at z > 2 have shown the rapid 
evolution of the line number density (Lu, Wolfe & Turnshek 
1991; Bechtold 1994; Kim et al. 1997; KCD). Weymann et al. 
(1998) from the HST QSO absorption line key project, how- 
ever, have shown that the redshift evolution of the line num- 
ber density is much more gradual at redshifts below z ~ 1.5 
than above. These results suggest that the transition be- 
tween the two different evolutionary rates occur somewhere 
in the range z ~ 1.3 — > 1.7. From numerical simulations, 
Theuns et al. (1998) and Dave et al. (1999) have demon- 
strated that the change in the evolutionary slope occurs at 
z ~ 1.7 due to the decreasing UV background at z < 2, 
assuming a QSO-dominated background. Recently, from the 
study of the Lya forest at 1.5 < z < 2.4, KCD conclude that 
the Lya forest at 1.5 < z < 4 shows a monotonic, continu- 
ous evolution with z, in terms of the line number density, the 
mean H I opacity and the correlation strength, both from the 
profile fitting analysis and from the optical depth analysis. 
KCD also show that the change in the number density evo- 
lution occurs at z ~ 1.2, suggesting that a QSO-dominated 
UV background used in numerical simulations underesti- 
mates the UV background and/or the enhanced structure 
formation at z < 2.5 (Theuns et al. 1998; Weymann et al. 
1998; Dave et al. 1999; Bianchi et al. 2001). Therefore, it is 
of importance to investigate the redshift range where dn/dz 
starts to change and to study any variations from sightline to 
sightline from larger statistical sample in order to constrain 
the results from numerical simulations and the nature of the 
ionizing sources. 

Here, we present the analysis of six QSOs covering the 
Lya forest at 1.5 < ZLya < 2.4 as well as two QSOs at 
3 < z Lya < 3.6 from high resolution (R ~ 45 000), high S/N 
(~ 35-50) spectra obtained with the VLT/UVES. This anal- 
ysis is an extension of the one by KCD with a larger uniform 
sample. We consider the evolutionary behavior of the Lyman 
forest systems over the redshift range 1.5 < z < 2.4, in par- 
ticular the line number density evolution and the clustering, 
and compare the results with those found by others at higher 
and lower redshifts. We have adopted the traditional Voigt 
profile fitting method to analyse the absorption lines in two 
ways. In the first approach, we have fitted the Lya forest 
region using the Lya lines only, since most previous studies 
have adopted this approach. In the second approach, we have 
fitted the Lya forest with the higher order lines of Lyman 
series (predominantly Ly/3) at wavelengths down to the ob- 
servational limit, 3050 A, for two QSOs at z om ~ 2.4 among 
8 QSOs presented in this study. This analysis includes the 
Lya forest from lower z in the regions where these higher 
order lines are present. In Section 2, we describe the UVES 
observations, data reduction and the Voigt profile fitting. In 



Section 3, we present the analysis of the Lya forest from the 
fitted line parameters, such as the number density. The Ly/3 
forest region is discussed in Section 4. The conclusions are 
summarized in Section 5. 



2 OBSERVATIONS AND DATA REDUCTIONS 

Table |l| lists the observation log for the QSOs observed 
with the VLT/UVES. In addition to HE0515-4414, J2233- 
606 and HE2217-2818 from KCD, we have included five 
more QSOs. Note that we extended the wavelength cover- 
age for the Lya forest towards J2233-606 and HE2217-2818 
in this study compared to those from KCD. The spectra 
were reduced with MIDAS/UVES and the resolution is about 
R ~ 45 000. The S/N varies across the spectrum and the typ- 
ical S/N in the Lya forest is ~ 40-50 for all the QSOs ex- 
cept for HE1347-2457, which is somewhat lower, ~ 35. The 
spectra were normalized locally with the 5th and 7th-order 
polynomial fitting. In order to avoid the proximity effect, we 
only include the Lya forest 3 000 km s _1 shortward the Lya 
emission in this study. See KCD for the details of the data 
reduction. 

Traditionally, the Lya forest has been thought of 
as originating in discrete clouds and has thus been an- 
alyzed as a collection of individual lines. These absorp- 
tion lines are generally fitted by Voigt profiles. From the 
Voigt profile fits, three line parameters (the absorption red- 
shift, z, the Hi column density, Nui in cm~ 2 , and the 
Doppler parameters, b in km s _1 ) are derived for each 
cloud. We have used the VPFIT program (Carswell et al.: 
http://www.ast.cam.ac.uk/~rfc/vpfit.html) to fit the lines, 
and for blended systems we have added the minimum num- 
ber of component clouds to ensure that the reduced x 2 is 
below an adopted threshold value of 1.3. Voigt profile fit- 
ting is not unique (cf. Kirkman & Tytler 1997; KCD), but 
we have fitted all the absorption lines in a consistent manner 
here to allow comparisons of the properties between different 
redshifts and sight-lines. 

We have adopted two approaches to define a sample of 
line parameters from VPFIT. Most previous studies have 
analysed the Lya forest longward the Ly/3 emission line to 
avoid the confusion with the Ly/3 forest from higher redshift 
absorbers. Thus, for comparison with other studies, we have 
used only the Lya absorption lines longward the Ly/3 emis- 
sion, in our first approach. All the QSOs in Table [i] have been 
fitted in this way, without considering their Ly/3 absorption 
profiles. The line parameters fitted with only the Lya pro- 
files define Sample A. We restrict our analysis to systems 
with /Vhi > 10 12 ' 5 cm -2 , the greatest value for the detection 
limit for all the quasars studied. At Nm = 10 12 ' 5 cm -2 , the b 
values in general span from 15 to 45 km s~ 2 where the larger 
value is effectively set by the detection limit. In addition, we 
have excluded the absorption lines with iVm > 10 17 cm -2 
in order to avoid Lyman limit systems (Note that KCD 
analyses the forest at Nm = 10 12,5-16 cm -2 . Since there 



are very few lines with /Vhi 



10 1 



cm , their results 



would be very similar even if the Nui range were changed 
to Nm = 10 12 ' 5 - 17 cm~ 2 ). 

Note that there is one high column density system with 
iVm = 10 1746 cm -2 at z = 3.192 towards Q0055-269 from 
the Lya absorption profile. There is, however, no corre- 
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Table 1. Observation log 



Aon 
QSO 


B a 


2em 


Wavelength 


Exp. time 


Observing Date 








( A ) 


(sec) 




HE0515-4414 


14.9 


1.719 


3050-3860 


19000 


Dec. 14, 18, 1999 


Ql 101-264 


16.0 


2.145 


3050-3870 


23400 


Feb. 10-16, 2000 


J2233-606 


17.5 


2.238 


3050-3860 


16200 


Oct. 8-12, 1999 








3770-4980 


12300 


Oct. 10-16, 1999 


HE1122-1648 


17.7 


2.400 


3050-3870 


26400 


Feb. 10-16, 2000 








3760-4975 


27000 


Feb. 10-16, 2000 


HE2217-2818 


16.0 


2.413 


3050-3860 


16200 


Oct. 5-6, 1999 








3288-4522 


10800 


Sep. 27-28, 1999 


HE1347-2457 


16.8 


2.617 


3760-4975 


18000 


Feb. 10-16, 2000 


Q0302-003 


18.4 


3.281 


4806-5771 


20000 


Oct. 12-16, 1999 


Q0055-269 


17.9 


3.655 


4634-5600 


18100 


Sep. 20-22, 2000 








4790-5729 


17000 


Sep. 20-22, 2000 



a Taken from the SIMBAD astronomical database. The magnitude of HE1347-2457 is from NED. 
Table 2. Analyzed QSOs 



QSO 


AA 


z LyQ dX a # of lines b 




<z> 


Sample A 



HE0515-4414 


3080-3270 


1.53- 


-1.69 


0.408 


63 


0-086%%9 B1 


2.1 


Q1101-264 c 


3230-3400 


1.66- 


-1.80 


0.381 


62 


0.085%%^ 






3500- 


-3778 


1.88- 


-2.08 


0.685 


99 


0.106%%^ 






3230-3400, 


3500-3778 


1.66- 


-2.08 


1.066 


161 


na 6 0.049 


2.1 


J2233-606 d 


3400- 


-3890 


1.80- 


-2.20 


1.209 


166 


n i c: fi 0.049 
U.10D_ Q51 


2.1 


HE1122-1648 


3500- 


-4091 


1.88- 


-2.37 


1.518 


234 


n 14fi 0.049 
U.14D_ Q51 


2.1 


HE2217-2818 


3510- 


-4100 


1.89- 


-2.37 


1.519 


214 


n ioqO.049 


2.1 


HE1347-2457 


3760-4335 


2.09- 


-2.57 


1.575 


233 


14 n0.049 


2.1 


Q0302-003 


4808- 


-5150 


2.96- 


-3.24 


1.152 


167 


n oo 4 0.049 


3.3 


Q0055-269 


4852- 


-5598 


2.99- 


-3.60 


2.638 


419 


o 401O.O49 
u -^ zl -0.051 


3.3 


Q0000-263 e 


5450-6100 


3.48- 


-4.02 


2.540 


312 


n 7-3Q0.049 
u. (oo_ 051 


3.8 



Sample B 



HE1122- 


-1648 


3200-3500 f 
3100-3500 


1.63-1.88 
1.55-1.88 


0.680 
0.893 


2is 


0.084^ 


HE2217- 


-2818 


3200-35 10 f 
3100-3510 


1.63-1.89 
1.63-1.89 


0.704 
0.917 


34S 


0.142%°^ 


The Ly/3 forest 


HE1122- 


-1648 


3674-4091 


2.02-2.37 


1.09 


195 




HE2217- 


-2818 


3674-4100 


2.02-2.37 


1.12 


190 





a For q = 0. 

b Over the column density range, Nm = 10 12,5-17 cm -2 . 

c There is a damped Lyo system at z = 1.8386. The Lya forest regions closer to the damped Lyet system by less than 50 A at each side 

have been excluded in the study. 

d See also Cristiani & D'Odorico (2000). 

e Lu ct al. (1996). 

f For the spectral region 3100—3200 A, the UVES spectra show a sharp decrease in S/N, so the continuum level is difficult to determine 
there. Consequently we exclude this region from the mean optical depth determination. Line counting applies for higher Hi column 
densities, where such uncertainties are less important, 
s For jV H , = 10 13 - 64 - 17 cm- 2 . 
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sponding Lyman limit, indicating that this high-A^i value 
is not real. In fact, the Ly/3 and the Ly7 absorption profiles 
corresponding to this system show a complex of at least 2 
lower- Ah i clouds, which is difficult to de-blend due to se- 
vere line blending from lower- z Lya forest. However, since 
we define Sample A from the results of fitting Lya absorp- 
tion profiles only, we treat this system as one with high Ah i 
and so have omitted it from the sample. 

Note also that for the optical depth analysis we have 
used all the available wavelength ranges between the Lya 
and Ly/3 emission lines without excluding any high- Ah i re- 
gions, except for Q1101-264. There is a damped Lya system 
at z — 1.8386 towards this QSO, and in this case we exclude 
the regions within 50 A of the central wavelength of this 
feature. 

In the another second approach, we have used the en- 
tire spectrum down to 3050A, fitting the Lya lines with the 
higher orders of Lyman series and adding the lower- z Lya 
lines simultaneously. The S/N at the Ly/3 forest is usually 
much lower than that at the Lya forest and the continuum 
fitting becomes more uncertain. It is of not useful to include 
the weaker Lya forest lines shortward the Ly/3 emission since 
its lower S/N would degrade the results from the statistical 
analysis. It is of importance, however, to include the cor- 
responding Ly/3 lines to determine JVm and b of saturated 
Lya fines at Ahi > 10 14 ' 5 cm -2 more reliably. Some of sat- 
urated lines are found to break into several lower- A^hi lines. 
Including these newly recognized lines of the Lya forest at 
Ahi > 10 14 ' 5 cm~ 2 is important if we are considering the ab- 
sorbers as individual entities for line counting and clustering 
analyses. 

During this process, it was found that sometimes the 
continuum needed to be re-adjusted by small amounts to 
obtain a satisfactory fit for all Lyman series including all 
the identified metal lines. Since fitting the Ly/3 forest region 
requires a higher S /N to determine a reliable continuum, not 
all the QSOs in Table |l| are suitable for this fitting method. 
In addition, the rapidly increasing numbers of the forest 
lines with z makes it difficult to fit all the spectral regions 
simultaneously at z > 3. For these reasons, we have selected 
two QSOs, HE1122-1648 and HE2217-2818, out of 8 QSOs 
in the sample. Both QSOs have S/N of 30-40 at 3200-3500 
A and of 10-30 at 3100-3200 A, suitable to explore the Ly/3 
forest at 2 < z < 2.4. 

In addition to the Ly/3 forest, this second fitting ap- 
proach provides the Lya forest at lower z, i.e. z < (1 + 
Zem) x (1025.72/1215.67) - 1. Since the Lya forest in the 
Ly/3 regions have a lower S/N, we have restricted this Lya 
forest in the Ly/3 forest region only for the study of higher 
Ahi forest at A?hi > 10 13 cm -2 , i.e. the line number density 
and the mean Hi opacity, in order to increase our statis- 
tics on z < 2.5. This Lya forest in the Ly/3 regions defines 
Sample B in this study. 

Keep in mind that all the analyses of the Lya forest here 
are from Sample A only, except the line number density and 
the mean Hi opacity. In any figures in subsequent sections 
filled circles represent Sample A, while filled squares are for 
Sample B. Table |^ lists the QSOs defined Sample A, Sample 
B and the Ly/3 forest. 

Metal lines were identified and removed as described in 
KCD. Different transitions of identified metal lines were fully 
taken into account when fitting the Lya forest. Metal lines 



could be assumed to be almost fully identified at z < 2.5, 
especially towards HE1122-1648 and HE2217-2818 where 
the analyses cover the full wavelength range down to 3050 

A. The only possible exception is for HE1347-2457, where 
incomplete coverage of the Lya forest may have resulted 
in some heavy element lines being missed. At higher-z, 
the identification of metal lines becomes more problematic 
due to severe line blending. Metal contaminations, however, 
should be less than 5 per cent at all redshifts. 

Figs. [jj, ^, ^, ^ and |E| show the Lya forest used in Sam- 
ple A: Q1101-264, HE1122-1648, HE1347-2457, Q0302-003 
and Q0055-269, respectively. The spectra are superposed 
with the fitted spectrum from the Voigt profile analysis (the 
sample A fitted line lists from the Voigt profile analysis 
with their errors are available electronically at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5). An ex- 
ample of the line lists is shown in Appendix in case of Q1101- 
264. The line lists of HE1 122-1648 and HE2217-2818, in- 
cluding the Ly/3 and Ly7 lines down to 3050A, i.e. Sample 

B, will be published elsewhere). The tick marks indicate the 
center of the lines fitted with VPFIT and the numbers above 
the bold tick marks indicate the number of the fitted line in 
the line lists. 



3 THE LYa FOREST 

3.1 The evolution of the line number density 

The line number density per unit redshift is defined as the 
number of the forest above a given A^hi per unit redshift. 
It is empirically described by dn/dz = (dn/dz)o{l + z) 1 , 
where (dn/dz)o is the local comoving number density of the 
forest. For a non-evolving forest in the standard Friedmann 
universe with the cosmological constant A = and the con- 
stant UV background, 7=1 and 0.5 for qo — and 0.5, 
respectively. Note that the measured 7 is dependent on the 
chosen column density thresholds, the redshift ranges and 
the spectral resolution (Kim et al. 1997; KCD). 

Fig. |^ shows the number density evolution of the Lya 
forest in the interval Ahi = ^q13.64-i7 cm -2 f rom g am _ 
pie A (filled circles) and Sample B (filled squares). This 
threshold has been chosen to be comparable to the equiv- 
alent width threshold of 0.24 A from the HST QSO ab- 
sorption line key project (Weymann et al. 1998), assuming 
AT H i = 1.33 X 10 20 W / /A§/, where W is the equivalent width 
in angstrom, Ao is the wavelength of Lya in angstrom, and 
/ is the oscillator strength of Lya. 

The dashed line represents the maximum-likelihood fit 
to the UVES and the HIRES data at z > 1.5. The line num- 
ber density of HE2217-2818 at 2 ~ 1.8 (Sample B) is higher 
than that of HE1 122-1648 at the same z range. On the other 
hand, the triangle at z ~ 1.55 from the HST observation of 
UM18 (which is considered as an outlier by Weymann et 
al. 1998) is well fit to the power-law derived at z > 1.5. As 
pointed out by KCD, (dn/dz) decreases as z decreases at 
1.5 < z < 4 with a consistent pattern and the change in 
dn/dz occurs at z < 1.5. 

The lower resolution (R ~ 1,000) of the HST/FOS 
data, however, makes it difficult to compare the results by 
Weymann et al. (1998) to the results from higher resolution 
ground-based observations. In addition, it should also be 
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Figure 1. The spectrum of Q1101-264 superposed with the fitted spectrum from the Voigt profile fitting. The residuals (the differences 
between the observed and the fitted flux) shown in the bottom part of each panel are shifted by —0.25. 
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Figure 2. The spectrum of HE1122-1648 superposed with the fitted spectrum from the Voigt profile fitting. The residuals (the differences 
between the observed and the fitted flux) shown in the bottom part of each panel are shifted by —0.25. 
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Figure 3. The spectrum of HE1347-2457 superposed with the fitted spectrum from the Voigt profile fitting. The residuals (the differences 
between the observed and the fitted flux) shown in the bottom part of each panel are shifted by —0.25. 
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Figure 4. The spectrum of Q0302— 003 superposed with the fitted spectrum from the Voigt profile fitting. The residuals (the differences 
between the observed and the fitted flux) shown in the bottom part of each panel are shifted by —0.25. 



noted that there is no reliable conversion from the equivalent 
width, W, to JVhi, without knowing the Doppler parameter, 
b, of the absorption lines. The conversion law between W 
and JVhi used in Fig. ^ is correct only if rm < 1. Due to 
the low resolution and to the large uncertainty in the con- 
tinuum fitting, many absorption lines from the HST QSO 
absorption line key project are unresolved and so Nhi may 
be underestimated. In the literature, b of 25 km s -1 is often 
assumed, under which the 0.24 A threshold corresponds to 
7Vhi = 10 14 cm -2 (cf. Savaglio et al. 1999; Penton, Shull & 
Stocke 2000; but see also Dave & Tripp 2001). 

Fig. ^ is the same as in Fig. [| except the A^hi range 
in dn/dz, Nm = 10 14-17 cm -2 . The dashed line is the 



Table 3. The line number density evolution 



range 


Median AT Hl 


(dn/dz)o 


7 


(cm -2 ) 


(cm -2 ) 






1Q13.1-14, a 


1Q 13.5 


50.0 


1.18 ±0.14 


1Q13.1-14, b 


10 13.5 


35.0 


1.42 ±0.16 


1q13. 64-17 


10 14.0 


6.1 


2.47 ±0.18 


10 14-17 


10 14.3 


1.8 


2.90 ±0.25 


10 14.5-17 


10 15.0 


0.5 


3.11 ±0.42 



a All the QSOs included. 

b Sample B of HE2217-2818 and HS1946T7658 are excluded. 
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Figure 5. The spectrum of Q0055-269 superposed with the fitted spectrum from the Voigt profile fitting. The residuals (the differences 
between the observed and the fitted flux) shown in the bottom part of each panel are shifted by —0.25. 
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Figure 5 - continued 



maximum-likelihood fit for the UVES and HIRES obser- 
vations. The points indicated by pentagons in Fig. ^ are 
taken from the Savaglio et al. (1999) analysis of the spec- 
trum of J2233-606 in the Ly/3 and the Ly-y regions. There 
are two important observational results evident from this 
figure. First, the Lya forest at 1.5 < z < 2 shows lower 
dn/dz than at z < 1.5, except the HE2217-2818 forest. Al- 
though the HE0515-4414 forest and the HE1122-1648 for- 
est could be considered to be consistent with the dn/dz at 
z ~ 1 within 2a, the Q1101-264 forest shows more than 3a 
difference. Since the universe expands and the overdensity 
evolution is rather smooth as a function of z for overdensi- 
ties corresponding to most Lya forest clouds, it is difficult 
to understand the higher dn/dz at z < 1.5 in terms of the 
overdensity evolution, in particular, a sharp transition at 
z ~ 1.5 shown in Fig. Instead, it is likely to be caused 
by the incorrect conversion between W and A^h i as well as 
the effect from the different resolutions. In reality, dn/dz 
should be somewhere in between Fig. ^ and Fig. |^ 

Second, the number density of the HE2217-2818 forest 
at z ~ 1.8 (Sample B) corresponds to those of UM18 and 



t Since the universe expands and the higher overdensity region 
collapses to form a structure, the same JVh i at different z samples 
different overdensity. The overdensity corresponds to Nn l ~ 



10 12 ' 7 cm~ 2 , N Kl ~ 10 13 ' 4 cm" 2 , and 7V H i 
2, 3, and 4, respectively (Schaye 2001). 



10 1 



at z 



J2233-606 at similar z. Even though the dn/dz of the UM18 
forest is discarded due to lower resolution, the HE2217-2818 
forest agrees well with the J2233-606 forest at z ~ 1.8. In 
fact, the higher dn/dz from HE2217-2818 and J2233-606 
is because that their sightlines include several higher- JVhi 
systems at iVm > 10 14 ' 5 cm -2 . This indicates that the lines 
of sight at z < 2.5 might not be as homogeneous as at z > 
2.5. 

In Fig. the dashed line is the maximum likelihood 
power law fit for z > 2.5 which has been extrapolated to 
lower redshifts. The dotted line gives a maximum likelihood 
fit for z < 2.5. A single power law does not give a satis- 
factory fit over the entire range (the dot-dashed line). Be- 
cause of the small number of sightlines and the different 
dz ranges for each sightline, it is not obvious how to inter- 
pret the dn/dz vs z results for this column density range. 
The general trend may be for dn/dz to decrease with de- 
creasing z, with inhomogeneity accounting for some of the 
higher dn/dz values at lower redshifts. Alternatively, dn/dz 
could be almost constant for z < 2.5 with a large scatter 
from sightline to sightline. Theuns et al. (1998) predict that 
dn/dz for the same column density range would show a non- 
evolution at z < 2. They do not, however, predict any spatial 
variation of dn/dz. 

When extrapolated from the behaviour of dn/dz at 
lower column density range, however, the former explanation 
gives a better interpretation of Fig. |^. In short, the higher 
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Ahi systems evolve more rapidly and show more evolved 
structures in their distributions along different lines of sight. 
Fig. ^|is the same as in Fig. |(| except the JVh i range used 



is Nhi = 10 



cm . The dashed line is the maximum- 



likelihood fit for the UVES and HIRES observations. As 
pointed out by KCD, the dn/dz at z ~ 2.7 (the diamond: 
HS 1946+7658) by Kirkman & Tytler (1997) is likely due 
to an over-fitting of higher S/N data than that of the rest 
of the observations. In fact, its thi fits within the fitted 
power- law (see Fig. [ll]) . On the other hand, the filled square 
at z ~ 1.8 more than 3a above the dn/dz at similar z is 
from the HE2217-2818 forest. It also shows a higher thi 
and a higher dn/dz than the fitted power-law due to several 
high Ahi forest. The dotted line represents the maximum- 
likelihood fit when Sample B of the HE2217-2818 forest and 
the HS1946+7658 forest are excluded. 

In Fig. ^, we assume that line blending at this Ahi 
range is not severe. Line blending, however, could lead to 
underestimate the line number density by as much as ~ 30 
per cent at z ~ 3 (see Hu et al. 1995; Giallongo et al. 1996). 
This would cause in part a slower evolution of line number 
density at z > 3. Although there is a line of sight with a 
different behavior, the lower Ahi forest also evolves with z 
monotonically at 1.5 < z < 4 in terms of the observed line 
density. 

Fig. [To] shows 7 as a function of the median Ah i for the 
different column density ranges used. Squares and diamonds 
represent the upper A^hi threshold as Ahi = 10 17 cm -2 and 
Nui = 10 16 cm -2 , respectively. Due to a lower number of 
the forest at JVhi = 10 16-17 cm -2 , 7 does not change signif- 
icantly with the upper JVhi threshold. The 7 values increase 
as the median Ahi increases (Lu et al. 1991; Weymann et al. 
1998). This result could be explained by two different sce- 
narios: blending and incompleteness of lower Ahi at higher 
z, and a change in the intrinsic properties of the Lya for- 
est. Keep in mind that the evolution of the UV background 
would predict the same 7 for all different column density 
ranges (Dave et al. 1999). 

Lower- Ah 1 lines are lost at higher z due to more severe 
line blending, i.e. incompleteness, resulting more line loss 
at higher z than at lower z. As a result, 7 could become 
smaller when lower- A^hi lines are included in the analysis 
(Giallongo et al. 1996; Kim et al. 1997; Dave et al. 1999). 
When available lines decrease at Aui.th > 10 14 ' 5 cm -2 , 7 

-,13.5-14.5 „„ 



10* 



7 



becomes very uncertain. At JVm.th 
shows a value of ~ 2.3. 

Unfortunately, quantifying the amount of line blending 
is not straightforward. The incompleteness corrections con- 
sidered in the literature have assumed a single power-law 
of the column density distribution. At lower- A^hi, there is a 
trend of higher power-law index as z decreases. This could 
not only due to incompleteness but also due to the real struc- 
ture evolution. At higher- A^hi at Ahi > 10 14 cm -2 , the on- 
set of rapid, non-linear collapse results in a deviation from a 
single power-law in the column density distribution, which 
is a function of z. Extensive Monte-Carlo simulations are 
required to establish the level of line blending as a function 
of z because there are large uncertainties in the continuum 
adjustments and varying S/N across the observed spectra. 
At the same time, these results should be compared with 
the ones from numerical simulations in order to constrain 
the significance of blending on the 7-median A^h i relation. 
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Figure 6. The number density evolution of the Lya forest over 
the column density range A^hi = 10 13,64— 17 cm -2 , which is com- 
parable to the HST data (open triangles) of Weymann et al. 
(1998). The data are shown for the binned sample for display. 
The filled symbols are derived from the UVES observations. Open 
squares, the star, open circles, and the diamond are taken from 
the HIRES data by Hu et al. (1995), Lu et al. (1996), Kim et 
al. (1997), and Kirkman & Tytler (1997), respectively. The hori- 
zontal error bars represent the z interval over which the number 
density was estimated. The vertical error bars represent the Pois- 
son lcr error. 



On the other hand, the increase of 7 as a function of the 
median Ah i could be a result of the intrinsic evolution of the 
Lya forest itself. Weaker Lya forest arising from lower den- 
sity gas expands faster than stronger Lya forest. Thus the 
fractional cross section to detect weaker Lya forest increases 
as z decreases (Dave et al. 1999). 

It is obvious, however, that a stronger variation in 
dn/dz for higher Ahi clouds also plays a part in the 7- 
median AThi relation (see Fig. |j). Due to a couple of higher 
dn/dz at z < 2.5 for higher column density forest, the ex- 
ponent from the maximum likelihood fit becomes smaller. 



3.2 The evolutionary behavior of dn/dz 



As shown in Figs. |6|-j9j, the break point in dn/dz may oc- 
cur at redshifts as low as z ~ 1, and very likely z < 1.5, 
rather than z ~ 1.7-2 as suggested by Theuns et al. (1998) 
and Dave et al. (1999). This shows that the UV back- 
ground does not decrease as rapidly as the QSO-dominated 
UV background by Haardt & Madau (1996) used in sim- 
ulations, and so galaxies could be a main source of the 
UV background at z < 2.5 and/or that the structure for- 
mation could be underestimated in simulations. Dave et 
al. (1999) show that dn/dz oc (1 + z) 5 ' 3 " 6 - 5 for a con- 
stant UV background in the flat Friedmann universe. For 
Am = 10 13 ' 64 " 17 cm" 2 at 1.5 < z < 4, /3 = 1.74T0.02. This 
indicates dn/dz oc (1 + z) 2 ' 25 , which is in agreement with the 
observations, dn/dz oc (1 + z) 2,44±0 18 . This result suggests 
a non-decreasing UV background at 1.5 < z < 4, unlike 
the decreasing one expected from the Haardt & Madau UV 
background at z < 2.5, i.e. the QSOs alone might not be 
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Figure 7. The number density evolution of the Lya forest over 
the column density range TVhi = 10 14— 17 cm -2 . All the symbols 
have the same meaning as in Fig. a 



Figure 9. The number density evolution of the Lya forest over 
the column density range Nn 1 = 10 131— 14 cm -2 . All the symbols 
unlisted in the figure are the same as in Fig. pi 
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Figure 8. The number density evolution of the Lya forest over 
the column density range Nn 1 = 10 14 ' 5- 17 cm -2 . All the sym- 
bols have the same meaning as in Fig. ^. The dot-dashed line 
represents the single maximum likelihoold fit over the entire red- 
shift ranges (see Table 0). The dotted line indicates the maximum 
likelihood fit for z < 2.5: dn/dz = 9f.3(l + 2 )-i-55±i.49 The 
dashed line represents the maximum likelihood fit for z > 2.5: 
dn/dz = 0.1 (1 -|-2:) 4 24±0 - 88 . The redshift 2.5 is chosen arbitrary. 



enough to explain the dn/dz evolution of the forest (KCD; 
Bianchi et al. 2001). 

Besides the break point in dn/dz at z ~ 1, Lu et al. 
(1991) suggest that dn/dz could be fit better with a double 
power-law at z ~ 2.3 for 1.6 < z < 3, using the minimum 
W threshold W = 0.36A. The minimum W threshold corre- 
sponds to JVhi = 10 14 7 , 10 14 4 and 10 14 2 cm" 2 for b = 25, 
30 and 35 km s" 1 , respectively. This threshold is similar to 
the one used in Fig. B. As clearly seen in Fig. @, there is no 



Figure 10. The value 7 as a function of the median JVhi 
at the different column density ranges. The upper Afjji thresh- 
old is fixed to be TVjji = 10 17 cm -2 (squares) and -/Vjji = 
10 16 cm -2 (diamonds), respectively. The lower iVni threshold for 
each median Af H 1, ^Hi.mcd = 10 13 3 , 10 13 6 , 10 13 9 , 10 14 3 , and 
10 15 cm -2 , corresponds to be 7V H i = 10 12 5 , 10 13 , 10 13 5 , 10 14 , 
and 10 14 ' 5 cm -2 , respectively. 

clear trend of a different evolution of dn/dz at higher and 
lower redshifts. 

The general trend shown in dn/dz confirms that dn/dz 
continues its evolution in the same general manner at least 
for 1.5 < z < 4, which indicates the importance of the Hub- 
ble expansion over the structure formation at this redshift 
ranges (Miralda-Escude et al. 1996; Theuns et al. 1998; Dave 
et al. 1999). Our new observations (Fig. ^|), however, show 
that the structure evolution becomes more important in 
some lines of sights at z < 2.5 and the clustering of stronger 
lines also increases as z decreases (Fig. [l6|) . 
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3.3 The mean Hi opacity 

As seen in Section 3.1, comparing the the line number den- 
sities at z > 1.5 with those at z < 1.5 is somewhat un- 
certain because of the different spectral resolutions used. 
The mean H I opacity provides more straightforward com- 
parisons between different qualities of data. Moreover, the 
mean H I opacity does not rely on the subjective line count- 
ing method, although it is more subject to continuum un- 
certainties (KCD). 

Fig. [ll] shows the thi measurements (the effective op- 
tical depth T e fr; exp _T ° ff = <exp _T >, where < > indicates 
the mean value averaged over wavelength), together with 
other opacity measurements compiled from the literature. 
Filled symbols are the thi measurements from the UVES 
data. The dotted line represents a widely used r c g from 
low-resolution data by Press, Rybicki & Schneider (1993). 
Table ^| lists the estimated r H i values from the UVES ob- 
servations (see KCD for the actual thi values at z > 1.5 in 
Fig. [l]} which is not listed in Table |). Table fi] lists the esti- 
mated thi values from HST observations, adding the equiv- 
alent widths. The different thi values of 3C273 at the dif- 
ferent z ranges by different studies suggest the uncertainty 
in deriving th i as well as a small-scale cosmic variance. 

The solid line represents the least-squares fit to the 
UVES (only filled circles and squares) and HIRES data: 
thi(z) = (0.0032±0.0009) (l+z) 3 ' 37±a20 . As stated in KCD, 
thi is well fit to a single power-law at 1.5 < z < 4 and the 
newly fit power-law is about a factor of 1.3 smaller than the 
Press et al. value from their lower resolution data. 

The square with higher thi at z ~ 1.8 is from HE2217- 
2818, which shows several high- Ah i clouds at Nm = 
10 14 ~ 17 cm -2 . This increases fm at z ~ 1.8, compared to 
that of HE1 122-1648 without high-A H i clouds at the same 
redshift range. The bold open circle at z ~ 2 is thi towards 
J2233-606, when two higher column density absorption sys- 
tems were excluded. Note that thi towards HE2217-2818 at 
z ~ 1.8 and J2233-606 are similar to the Press et al. value 
at the same redshifts. This confirms that the higher Press 
et al. value is in part due to the inclusion of high column 
density systems in their low resolution sample as well as 
the uncertainty in the continuum displacement, as noted by 
KCD. Keep in mind that most QSOs used in our analysis 
both from the UVES and Keck observations do show few 
high column density systems at Nai > 10 17 cm -2 . In the 
case of having a damped system along the sightline, such 
as Q1101-264 and Q0000-263, we excluded the regions in 
which the damped systems are located. 

There is a large uncertainty in thi from HST obser- 
vations due to their lower S/N, lower resolution, and lower 
absorption line densities with the presence of weak emis- 
sion lines, all resulting in an unreliable local continuum fit. 
This, in general, leads to a tendency to miss weak lines and 
so an underestimate thi- There is, however, indication of 
slow-down in the evolution of thi. Despite a large scatter 
at z ~ 0.1, the median thi at z ~ 0.1 is 0.019, a factor 
of 4.3 larger than the value extrapolated from at z > 1.5. 
Although open triangles (Weymann et al. 1998) might un- 
derestimate thi significantly, the different evolution in thi 
occurs at z < 1, a bit lower than that suggested from dn/dz. 

Dave et al. (1999) simulated the evolution of fm, as- 
suming different cosmologies and normalizing their results 




1 +2 



Figure 11. The Hi opacity as a function of z. Filled symbols 
represent the mean H I opacity from the UVES data. The open 
circle at <z> = 2.0 represents thi of J2233-606 when two high 
column density systems are excluded. Other symbols at z > 1.5 
are from: open circles (Hu et al. 1995), large square (Lu et al. 
1996), the open triangle (Kirkman & Tytler 1997), and diamonds 
(Rauch et al. 1997). Symbols at z < 1.5 with arrows are from: 
open diamonds (Impey et al. 1996), open triangles (Weymann et 
al. 1998), open stars (Impey et al. 1999) and open squares (Pen- 
ton et al. 2000). Due to the low resolution and low S/N spectra 
of the HST observations, r could be highly underestimated if the 
Lya forest at 7Vhi < 10 14 cm -2 contains the bulk of the neu- 
tral hydrogen, thi estimates, while the y-axis error bars were 
estimated from simply changing the adopted continuum by ±5%. 
The dotted line represents the commonly used formula by Press 
et al. (1993), t H i(z) = 0.0037(1 + z) 3A6 . The shaded area en- 
closed with dot-dashed lines indicates the ranges of f expected 
from different cosmological simulations by Dave et al. (1999). 

to observations at z ~ 3. The shaded area in Fig. ^ shows 
the ranges of thi from their simulations at z < 2. Although 
the estimation of thi from HST observatons could be un- 
derestimated by a large factor, it is not in good agreement 
with the simulated thi. In particular, at z ~ 0.5, there are 
some observations which show lower th i than the simulated 
one by more than a factor of 10. This could indicate that 
the QSO-dominated UV background by Haardt & Madau 
(1996) used in the simulations does not represent the real 
UV background at z < 2, i.e. this Haardt & Madau QSO- 
dominated background underestimates the real UV back- 
ground at z < 2, as suggested by the evolution of dn/dz. 
In addition, the structure in the Lya forest becomes more 
highly patchy at z < 2 than the simulations predict (see 
Fig. |8J). In the latter case, some lines of sight could have a 
very lower rm. 

3.4 The differential density distribution function 

The differential density distribution function, /(JVhi), is de- 
fined as the number of the absorption lines per unit absorp- 
tion distance path and per unit column density as a func- 
tion of Ah i. The absorption distance path X(z) is defined by 
X(z) = i[(l+z) 2 -l] forg = OorbyA(z) = §[(l+z) 3/2 -l] 
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Table 4. The mean H I opacity compiled from the literature 
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a The equivalent widths listed in the referenced paper are only 
for the core of the profile. Impey et al. (1996) gave a scale factor 
of 1.33 for isolated lines and of 1.18 for blended lines, to recover 
the total equivalent widths. We have used a scale factor of 1.255, 
a mean value. 

Ref.: 1. Bahcall et al. (1993); 2. Impey et al. (1996); 3. Impey et 
al. (1999); 4. Penton et al. (2000). 



for q — 0.5 in the standard Friedmann universe (see Table ^ 
for q — 0). Empirically, /(Am) is fitted to a power law: 
/(AT H1 ) = AiV-f . 

Fig. |l2| shows the observed /(Am) at different redshift 
ranges without the incompleteness correction due to line 
blending. The dotted line represents the incompleteness- 
corrected /(Am) at < z > = 2.85 from Hu et al. (1995), 
/(Am) = 4.9 x 10 7 TV", 1 ' 46 . Triangles for the damped Lya 
systems at 1.5 < z < 2.5 are taken from Storrie-Lombardi & 
Wolfe (2000). At z ~ 2.1, /(Am) at Am = 10 125 " 14 ' 5 cm" 2 
is in good agreement with the incompleteness-corrected 
/(Am) at z ~ 2.8. It is also true for /(Am) at z ~ 
3.3, although the goodness-of-the-fit is lower at Am = 
10 12,5 ~ 13 cm -2 . In general, /(Am) is well approximated to a 
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Figure 12. The differential density distribution functions at < 
2> = 2.1, 3.3 and 3.8, without the incompleteness correction. The 
data are shown in the binned sample for display. Open triangles 
are the differential density distribution functions for the damped 
Lya systems at 1.5 < z < 2.5 (Storrie-Lombardi & Wolfe 2000). 
The dotted line represents the incompleteness-corrected /(JVhi) 
at <z> = 2.85 from Hu et al. (1995), /(JV H i) = 4.9 X 10 7 TV" 1 ' 46 . 



single power-law f(Nm) oc iV^, 1,5 at 7V H i = 10 1 
(Petitjean et al. 1993). 

As noted by Petitjean et al. (1993) and Kim et al. 
(1997), /(JVhi) starts to deviate from the empirical power- 
law at Ah 1 > 10 14 cm" 2 . The amount of this deviation 
increases as z decreases since the higher Am forest disap- 
pears more rapidly as 2 decreases. In addition, the devia- 
tion Ah 1 at which /(Am) starts to deviate decreases as z 
decreases. At <z> = 3.8, 3.3 and 2.1, /(Am) deviates from 
the power-law at Am ~ 10 16 cm" 2 , Am ~ 10 14,5 cm" 2 and 
A H 1 ~ 10 14 - 2 cm" 2 , respectively. At Am > 10 15 ' 6 cm" 2 , 
the deviation from the single power-law increases more than 
3<r at z ~ 2.1 and more than 2er at z ~ 3.3. Table ^ lists 
log A and (3 for the various column density ranges from the 
maximum-likelihood fit. 

Fig. |l3| shows j3 as a function of the Am range. The 
x-axis error bars represent the Am range used in the fit. 
For Ah 1 < 10 14 ' 5 cm" 2 , (3 increases as Am increases, in 
part due to incompleteness at Am < 10 13 cm" 2 (Hu et al. 
1995; Giallongo et al. 1996, see also Section 3.1). For Am > 
10 14 ' 5 cm" 2 , (3 is ill-defined since there are not enough lines 
in the samples to get a reliable /3, especially at z ~ 2.1. 

Although it is statistically uncertain, there is a sugges- 
tion of increasing j3 as z decreases for a given Ah i range in 
Fig. [l3| Even if we exclude the lower- Ah 1 range distorted 
by line blending, the similar trends hold. Table |s| lists the 
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Figure 13. The slope /3 as a function of the JVhi range used in 
the maximum likelihood fit. 



P values at < 2 > = 0.03 using the conversion from W by 
Penton et al. (2000) and the /3 value at <z> =0.17 using 
the direct profile fitting analysis by Dave & Tripp (2001) 
for different column density ranges. Even if we discard the 
P values by Penton et al. (2000) due to their W analysis, 
the P value at < z > = 0.17 suggests that /3 increases as z 
decreases. This is another way of showing that the strong 
Lycv forest evolves rapidly at < z < 4. 

3.5 The differential mass density distribution of 
the Lya forest 

Schaye (2001) has shown that the shape of the differential 
density distribution function reflects that of the differential 
mass density distribution function as a function of the gas 
density. This differential mass density distribution function 
of the Lya forest is the mass density per unit overdensity 8. 
Assuming that the gas is isothermal, it can be described by 

dfL 



dlog(l + <5) 



7.6 xlO- 9 /*- 1 ^/ 3 ^ - 59 

f \ V3 
Jg 



N*{ 3 f(Nm,z), 



(1) 



where Ahi is related to 8 through 



A H , ~ 2.7 x 10 13 (1 + 5) 1 - 5 T 4 ~ ' 26 



/i^£\ 9 / 2 f n b h 2 

\ 4 J \ 0.02 



12 

3/2 



0.16, 



(2) 



(Schaye 2001). The parameter h is the Hubble constant di- 
vided by 100, the H I photoionization rate T = Ti 2 x 10 12 
s _1 , the temperature of the Lycv forest T = Ti x 10 4 K, 
/(Ah 1,2) is the differential density distribution function as 
a function of z, f g is a fraction of mass in the Lycv forest 
and fib is the baryon density (Schaye 2001). As the gas is 
assumed to be isothermal in Eq. 1, T is assumed to be con- 
stant (a = in T oc (1 + 8) a ), i.e. not a function of (1 + 8). 
The df2g/dlog(l+5) values, however, do not vary much with 
a. 



2001). If the differential density distribution function is a 
single power-law, then df2 g /dlog(l + 8) oc A^ 3- ' 3 . There- 
fore, the shape of dQ g /d\og(l + S) reflects the deviation from 
the single power-law of f(Nm, z). As seen in Fig. |l3[ P is not 
a constant. If P < 4/3, df2 g /dlog(l + <5) increases. Therefore, 
dfl g /d log (1 + 5) is expected to increase at smaller Ahi and 
to decrease at larger Nm- 

Fig. |l4| shows the differential mass density distribution 
function dQ g /d log(l + 8) as a function of Ahi at three red- 
shifts. The arrow in Fig. [l4| indicates the direction towards 
which dQ g /dlog(l + 6) moves if T4 or T12 increase, i.e. pre- 
serves the shape of dQ, g /dlog(l + S). The QSO-dominated 
UV background by Haardt & Madau (1996), Jhm, was as- 
sumed in the Friedmann universe with h = 0.65. Other pa- 
rameters are obtained from the observations, while T4 was 
read from figure 3 by Schaye et al. (2000). 

There is a turnover at Nm ~ 10 13 cm -2 due to line 
blending as seen in Fig. [12] and in part due to the real defi- 
ciency of these weaker lines (cf. Schaye 2001). The amount of 
the turnover is more significant at z ~ 3.8 and the turnover 
occurs at lower JVhi at lower z. At Nm ~ 10 13 5 " 15 cm -2 , 
df2 g /dlog(l + 8) is roughly approximated by a single power- 
law, with a slope being slightly steeper at lower z. This indi- 
cates that P is larger at lower 2 and that P at all z is larger 
than 4/3. At 7V H i > 10 15 cm -2 , df2 g /dlog(l + 6) is not 
clearly defined due to a larger size of the Nn 1 bin (also due 
to a smaller number of lines at this column density regime) . 

The three curves for df2 g /dlog(l + S) at different red- 
shifts agree with each other reasonably well with each other 
if one applies a redshift-dependent offset in Nm, at least in 
the range A^hi ~ 10 13-15 cm" 2 (referenced to < z >= 3.3). 
In fact, df2 g /dlog(l + 8) as a function of (1 + 8) does 
show a similar shape at different z. Eq. |i| shows that any 
change in z and Ti2 (thus TV) reflects the relation between 
Ah 1 and (1 + 8). Thus, Fig. indicates that the redshift- 
cvolution of the differential density distribution function 
seen at Am = 10 13-15 cm -2 reflects the change in V at 
different 2 (Schaye 2001). 



3.6 The two-point function of the flux 

The line width of the absorption lines (the b parameter from 
the Voigt profile fitting) provides valuable information on 
the thermal temperature of the absorption clouds, if the 
Lya forest is thermally broadened. In particular, the lower 
envelope of the Ah 1-6 diagram is interpreted to give an up- 
per limit on the temperature (through b) of the forest as a 
function of Ah i (Schaye et al. 1999; McDonald et al. 2000; 
Ricotti, Gnedin & Shull 2000; Schaye et al. 2000; KCD; Kim, 
Cristiani & D'Odorico 2002). However, it is not straightfor- 
ward to determine this minimum b value as a function of 
Ah 1, particularly because of profile fitting errors in blended 
features. 

Instead of using the fitted parameters, the direct use 
of normalized flux, F, has been proposed (Miralda-Escude 
et al. 1997; Bryan et al. 1999; Machacek et al. 2000; The- 
uns et al. 2000). Among these flux-based properties of the 
Lya forest, the two-point function of the flux provides the 
profile shape of absorption lines at 30 < Av < 100 km s -1 
(Machacek et al. 2000; Theuns et al. 2000). 



© 2002 RAS, MNRAS 000, 000-000 



16 T.-S. Kim et at 



Table 5. The power-law fit of the distribution functions, /(-/Vhi) = ANj 
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a Penton et al. (2000). 
b Dave & Tripp (2001). 
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Figure 14. The differential mass density distribution of the Ly. 
forest as a function of JVjj t . Among the parameters used in Eq. 
and Eq. Ti2 was read from figure 8 of Haardt & Madau (1996) 
and T4 was read from figure 3 of Schaye et al. (2000) for their 
simulations under Jjjm at <5 = 1. The parameters / g , f2 b /i 2 and 
h are assumed to be 0.16, 0.02 and 0.65, respectively, in the flat 
Friedmann universe. At < z > = 2.1, 3.3 and 3.8, (74,ri2) is 
(1.41,1.35), (1.57,1.07) and (1.62,0.7), respectively. The arrow 
indicates the direction towards which df2g/dlog(l + <5) moves if 
T4 or Ti2 increases. 



The two-point function of the flux, P(Fi,F 2 , Av), is the 
probability of two pixels with the Av velocity separation 
having normalized fluxes F\ and F%. It is described as 



AF(Av,5F!) 



5Fi 



(F 1 -F 2 )P(F 1 ,F 2 ,Av)dF 2 



dF 1 /SF 1 , (3) 



where AF(Av, 5Fi) is the mean flux difference between two 
pixels with F\ and F 2 , which are separated by Av (Theuns 

et al. 2000). 

Fig. |l| shows AF for 0.05 < F < 0.15 at different red- 
shifts. We measured the A«(0.3) value, following the defini- 
tion by Machacek et al. (2000), the width of AF(Av,6Fi) 
at which AF(Av, 8F1) becomes 0.3. The individual QSOs 
show a different Au(0.3) even at the same z. The measured 
Av(0.3) averaged at each 2 is 30.2 km s _1 , 33.7 km s _1 and 
48.8 km s _1 at < z > = 2.1, 3.3 and 3.8, respectively. This 
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Figure 15. The AF(Aw,c5Fi) profile for 0.05 < F < 0.15 at 
different redshifts. The left panel is for individual QSOs of Sample 
A, while the right panel is averaged for each z bin. 



result might indicate that the line profile becomes broader 
as z increases. In fact, Theuns et al. (2000) note that at a 
given z a simulation with a hotter gas temperature shows a 
wider AF(Av,8F\) profile than a simulation with a lower 
gas temperature. 

Decreasing A«(0.3) as decreasing z, however, is not an 
indicative of decreasing in temperature. The reason of de- 
creasing Av(0.3) is that AF(Av, 5Fi) must be asymptotic to 
F — 0.1 and F decreases as z increases (Theuns et al. 2000). 
Since F depends on z, comparing Aw(0.3) at different z di- 
rectly does not provide any information on the temperature. 

There is another way to look at decreasing Aw(0.3) with 
decreasing z. From the minimum cutoff b distribution at each 
z from the same data, Kim et al. (2002) show that the min- 
imum cutoff b values increase as z decreases (cf. Kim et al. 
1997). This result shows that the b (i.e. line width) values of 
absorption lines increases as z decreases. This result is con- 
trary to the result from AF(Av, SFi). This analysis of the 
minimum cutoff b values suggests that AF(Av, 5Fi) actually 
measures the relative amounts of line blending and higher- 
./Vhi forest, when it is compared at different redshifts. The 
Aii(0.3) values increase as z increases due to more severe 
blending and more higher- A^hi forest at higher z, not due to 
increasing temperature of absorbing clouds. AF(Av, 6F1) 
for higher F ranges, i.e. excluding higher column density 
lines, does not show any z-dependence. 
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3.7 The optical depth correlation function 

We analysed the clustering properties of the Lya forest, us- 
ing the step optical depth correlation function (Cen et al. 
1998; KCD). The optical depth correlation function is less 
biased than the line correlation functions since the fitted 
lines from Voigt profile fitting are not unique. 

The step optical depth correlation function £ T , S is de- 
fined by 

f T , (A .) = <r s (v + Av)r s (v)> ^ (4) 

\ T s S> 

where the step optical depth, r s , is r s = if r b s < Tmm and 

7"s — 1 if Tpb s ^ Tmin • 

Fig. |l6] shows £ T . s (An) averaged at each z for r m j n = 2. 
The Tmin = 2 corresponds to ATm ~ 10 13 ' 6 cm -2 , if b is 
assumed to be b = 30 km s _1 . Dotted lines, dashed lines 
and the dot-dashed line represent the individual QSOs in 
each z bin at < z > = 2.1, 3.3 and 3.8. The large error 
bars at < z > = 2.1 are due to the fact that the correla- 
tion strengths of different QSOs at similar redshifts show a 
wide range. Although there is a large scatter at <z> = 2.1, 
the step optical depth correlation functions show a strong 
clustering at v < 100 km s _1 . In addition, the step opti- 
cal depth correlation strength increases as z decreases in 
general. At < z > = 2.1, £ T , S (50 km s^ 1 ) shows a ~ 10<r 
significance compared to that at < z > =3.3. This stronger 
clustering of the Lya forest at z < 2 compared to higher z 
is expected from the differences in dn/dz and thi along the 
different sightlines. Studies of the two-point velocity corre- 
lation strength using the fitted line parameters at z < 1.5 
and z > 2.5 also show a velocity correlation at Av w 50-500 
km s~\ although its significance is much smaller than that 
of the step optical depth correlation function (Ulmer 1996; 
Cristiani et al. 1997; Penton et al. 2000). 



4 THE LY/3 FOREST 

At A?hi > 10 14 cm -2 , the Lya forest absorption lines start 
to saturate. For single lines, this saturation results in an 
increasingly inaccurate estimation of Ah i and b as the H I 
column density increases, at least up to the stage where 
damping wings become significant. More importantly, satu- 
ration also results in the uncertainty in the number of sub- 
components present. In reality, some high- Ah i clouds are 
not a really high- Ah i cloud, but a blend of several lower- 
A^hi clouds. Many of these saturated lines could have more 
accurate measurements of Ahi, b and a number of subcom- 
ponents if the higher order Lyman lines are fitted simulta- 
neously with the Lya, since the oscillator strengths decrease 
monotonically as one goes up the Lyman series. 

Fitting the absorption profiles with higher-order series 
becomes more time-consuming as z increases. In some cases, 
severe contamination by lower- z Lya forest and by the other 
higher-order Lyman series from higher- z forest make it very 
difficult to deblend the saturated Lya lines properly. In addi- 
tion, S/N of > 15-20 is required to determine the continuum 
with adequate reliability. For these reasons, we have only fit- 
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Figure 16. The step optical depth correlation functions at < 
z > = 2.1 (squares), < z > = 3.3 (diamonds) and < z > = 3.8 
(triangles). Error bars are the ltr statistical errors from the sample 
size. 

ted two QSO spectra with S/N ~ 10-30 in the Ly/3 regions, 
HE1 122-1648 and HE2217-2818§. 

4.1 The line number density evolution of the Ly/3 
forest 

Fig. |l^ shows the line number density, dn/dz, of the Ly/3 
forest (filled symbols) and the Lya forest both from Sample 
A and Sample B (open symbols). For AThi = 10 14 ~ 17 cm -2 , 
dn/dz of the Ly/3 forest is within the error bars of dn/dz of 
the Lya forest. In the case of HE2217-2818, dn/dz of the 
Ly/3 forest is less than lcr different from that of the Lya 
forest (14 versus 15 lines). For HE1122-1648, it becomes 27 
versus 25 lines, well within the error bars. On the other hand, 
for A?hi = 10 13,1 ~ 14 cm~ 2 , dn/dz of the Ly/3 forest systems 
shows a larger deviation. For HE2217-2818 (for HE1122- 
1648), it is 74 versus 68 lines (77 versus 62 lines). The main 
reasons for this difference in dn/dz come in part from the 
continuum re-adjustment in the Lya-Ly/3 fit (which has a 
more significant effect on the lower AThi forest) and in part 
from the more correct deblending of the higher- A%i forest. 
When the higher- ATm forest is deblended into more compo- 
nents, the number of lower- A^hi components in the forest 
increases. This results in the higher number of smaller A^hi 



§ Note that these two spectra have S/N of 20-30 in 3200-3500 A. 
Shorter than 3200 A, S/N decreases very rapidly, S/N ~ 10-20. 
We are, however, mainly interested in constraining a reliable ./Vh i 
determination of saturated ./Vh i clouds in the Lya forest regions 
and the higher- %, clouds in the Ly/3 forest regions. Therefore, 
lower S/N at 3100-3200 A does not limit to study the Ly/3 forest 
and the higher- Whi Lya forest in the Ly/3 regions. 
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Table 6. The power-law fit of the distribution functions, 
f(N Bl ) = AN-f at <z> = 2.2 
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Figure 17. The line number density evolution of the Lya forest 
and the Ly/3 forest. Open circles represent all the different sym- 
bols from Fig. ^, while open triangles are from Fig. ^j. The dashed 
line and the dotted line represent the maximum-likelihood fit of 
the Lya forest for each TVjji range in the fit. 



forest from using the line parameters determined from fit- 
ting the Lya- Ly/3 forest. 

In general, however, the uncertainty in dn/dz of the 
Ly/3 forest is in the same amount or less than the uncer- 
tainty resulted from the different sightlines, i.e. the cosmic 
variance. The determination of dn/dz from the Lya forest 
alone does not affect the overall results of dn/dz at least 
at 1.5 < z < 2.5. At higher z, more severe line blending is 
likely to produce higher uncertainties in the reliable number 
of subcomponents in the saturated lines. This might result 
in a significantly different dn/dz of the Ly/3 forest at z > 3, 
i.e. a higher dn/dz, compared to that of the Lya forest. 



4.2 The differential density distribution function 
of the Ly/3 forest 

Fig. [TJj shows the differential column density distribution 
function, /(TVhi), for the Lya-only fit (open squares) and for 
the Lya-Ly/3 fit (filled circles). At Nm < 10 143 cm" 2 , there 
is no significant difference in the Ly/3 forest and the Lya for- 
est. The Kolmogorov-Smirnov test gives the KS statistic of 
0.06 and the probability of 0.56. Note that a slight difference 
in f(Nni) is due to the continuum fitting adjustment when 
the Lya and the Ly/3 forest lines were fitted simultaneously. 

The difference seems to become more noticeable at 
Nm > 10 14 ' 3 cm -2 . This is, however, in part caused by 
the fact that there are not many high- Ah i clouds in the 
sample (22 lines in the Ly/3 forest sample and 18 lines in 
the Lya forest sample) and in part by the fact that Fig. |l8| 
shows the binned data for the display purpose. In fact, the 
Kolmogorov-Smirnov test gives the KS statistic of 0.22 and 
the probability of 0.65. Although it is a statistically small 
sample, Table H lists f{Nm) at different Ahi ranges. Within 
the error bars, /(Ahi) from both the Lya forest and the Ly/3 
forest is in good agreement. 
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7.40 ± 1.31 
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1.74 ±0.17 


Ly/3 


7.52 ± 1.34 


1.43 ±0.04 


11.92 ± 0.75 
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4.3 The distribution of the Doppler parameters of 
the Ly/3 forest 

One of the important physical properties of the Lya forest, 
the temperature of the absorbing gas, can be derived from 
the distribution of the Doppler (6) parameters. At higher z, 
the absorption lines are broadened by the thermal motion as 
well as the bulk motions. Therefore, the lower cutoff enve- 
lope of the Nui-b distribution constrains an upper limit on 
the thermal temperature of the forest as a function of Ahi 
at a given z (Hui & Gnedin 1997; Schaye et al. 1999; Mc- 
Donald et al. 2000; Ricotti et al. 2000; KCD; Dave & Tripp 
2001; Kim et al. 2002). 

Fig. jHj] shows the distribution of the Doppler param- 
eters of the Lya (dashed lines) and the Ly/3 forest (solid 
lines). For Nm = 10 12 - 5 " 14 ' 5 cm" 2 (the left hand panel), 
the median b is 27.6 km s" 1 (the Ly/3 forest) and 26.84 km 
s _1 (the Lya forest). There is no significant difference in the 
b distribution between the Lya forest and the Ly/3 forest for 
this Ah i range. This result confirms the validity of the de- 
termination of the lower cutoff b values using the Lya forest 
in previous studies since the lower cutoff has been derived 
at Nm = IO 12,5-14 5 cm -2 . Since these lines are not satu- 
rated, lower cutoff b values from the Lya forest would not be 
changed, when the Lya forest is fit simultaneously with the 
Ly/3 forest. In fact, the Nm~b diagram from the Lya forest 
and the Ly/3 forest does not show any significant difference. 

The right hand panel of Fig. FJ| shows the b distributions 
for Nm = 10 14 5 " 17 cm -2 . The median b is 34.0 km s" 1 
(the Ly/3 forest; 15 lines) and 32.7 km s" 1 (the Lya forest; 
13 lines). Although the median values are not statistically 
robust due to the small number of available lines, there is 
no significant difference between the Lya and the Ly/3 forest 
for this higher /Vhi range, either. 

Table M lists the observed b values and the median b val- 
ues from the Ly/3 forest as well as the ones at z ~ 0.15. Our 
characteristic b- values at z ~ 2.2 are very similar to those re- 
ported by Shull et al. (2000) from the curve-of-growth analy- 
sis. On the other hand, our values are somewhat larger than 
those of Dave & Tripp (2001) obtained from a profile fitting 
analysis. Table ^ shows that there is no significant differ- 
ence of b and median b values between the Lya forest and 
the Ly/3 forest at z ~ 2.2. Although b and median 6 values 
are dependent on Ah i ranges, Table |^ suggests that the me- 
dian 6 value does not increase, i.e. constant or decreasing, 
at z < 1.5 from at z ~ 2.1 for two Nm ranges considered in 
Table | (see KCD; Kim et al. 2002). 
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Table 7. The characteristic b values of the Ly/3 forest 



<z> 


(cm- 2 ) 


^median 

(km s- 1 ) 


b 

(km s- 1 ) 


Ref. 


0.15 


l 13 .76-17, a 


28 


31.4 ±7.4 


1 


2.2 


^OlS-76-17 


28.5 


30.2 


2 


0.17 


10 13-17 


22 


25 


3 


2.2 


10 13-17 


27.0 


32.6 


2 



a The equivalent width threshold W > 0.2A has been translated 

to JV H i = 10 13 - 76 cm- 2 , assuming b = 30 km s _1 . 

Ref.: 1. Shull et al. (2000); 2. This study; 3. Dave & Tripp (2001). 



Table 8. The characteristic b values of the Lya forest 



<z> 



(cm- 2 ) 



"median 

km s -1 



b 

km s _1 



Ref. 



0.15 


Iq13. 76-17, a 


28 


31.4 ±7.4 


1 


2.1 


Iq13. 76-17 


29.5 


32.7 


2 


3.3 


1Q13. 76-17 


27.5 


30.7 


2 


3.8 


Iq13. 76-17 


31.0 


39.5 


3 


0.17 


10 13-17 


22 


25 


4 


2.1 


10 13-17 


27.0 


31.5 


2 


3.3 


10 13-17 


28.1 


31.6 


2 


3.8 


10 13-17 


29.2 


35.9 


3 



a The equivalent width threshold W > 0.2A has been translated 

to Nn 1 = 10 13 - 76 cm -2 , assuming b = 30 km s — 1 . 

Ref.: 1. Shull et al. (2000); 2. This study; 3. Lu ct al. (1996); 4. 



Dave & Tripp (2001). 



5 CONCLUSIONS 

We have analyzed the properties of the Lya forest at 
Nm = 10 12,5 ~ 17 cm- 2 toward 8 QSOs, using high resolu- 
tion (R ~ 45 000), high S/N (~ 25-40) VLT/UVES data. 
The analyses presented here extend the ones by KCD, using 
the data obtained and treated uniformly. In addition, the 6 
QSOs at 1.5 < 2 < 2.5 enable us to explore the cosmic vari- 
ances along different sightlines. Combined with other high- 
resolution observations from the literature as well as HST 
observations at z < 1.5, we have studied the properties of 
the Lya forest as a function of z. 

To analyse the properties of the Lya forest, we have 
used two profile fitting approaches. In the first analysis, wc 
have only fitted the Lya absorption profiles in between the 
Lya and the Ly/3 emission lines. In the second analysis, we 
have fitted the Lya forest simultaneously with the higher or- 
der lines of the Lyman series down to 3050 A, for two QSOs 
at 2 cm ~ 2.4. The second analysis has been adopted to probe 
the properties of the Ly/3 forest at z ~ 2.2 and to investigate 
the general differences between studies from the traditional 
Lya-only fits and from the higher orders of the Lyman se- 
ries fits. In addition, we have also applied the optical depth 



analysis. For the lines with TVhi = 10 



cm , we have 



in general confirmed the conclusions by KCD derived from 
a smaller sample than that of this study. We have found: 

1) The line number density of the Lya forest, dn/dz, 
is fit well by a single power-law and shows a steeper evo- 
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Figure 18. The differential column density distribution function 
from the Lya-only fit and the Lya and Ly/3 fit. The data are 
shown in the binned sample for the display purpose. The dotted 
line represents the incompleteness-corrected /(-/Vjii) at <z> = 
2.85 from Hu ct al. (1995), /(%,) = 4.9 X 10 7 W Hl L46 . 




20 40 60 80 20 40 60 80 100 
b (km s" 1 ) 6 (km s" 1 ) 

Figure 19. The distribution of the Doppler parameters for 
the Lya and Ly/3 forest. The left hand panel is for JVhi = 
:K)12.5-14.5 cm -2 , while the right hand panel is for JV H , = 
1Q14 5-17 cm -2 gQiid lines are for the Ly/3 forest and dashed 
lines are for the Lya forest. 



lution at higher- JVhi forest at 1.5 < z < 4. For /Vhi = 

10 13.64-17 cm -2 ; dn/dz K (! + Z ) 2 -47±0.W p or jy H j = 

10 14 ~ 17 cm" 2 , dn/dz oc (1 + 2 ) 2 - 90±0 - 25 . The change in the 
evolutionary behavior of dn/dz occurs at z ~ 1. In addition, 
forest regions with under-abundant line density and regions 
with over-abundant line density start to appear at z < 2.5. 
This small-scale variation is more significant for higher /Vhi 
forest. This is probably caused by the small-scale structure 
evolution in the forest, i.e. the increasing clustering of high 
A^hi forest at z < 2.5. 



© 2002 RAS, MNRAS 000, 000-000 



20 T.-S. Kim et al. 



2) The mean H I opacity, t, is also well approximated by 
a single power law at 1.5 < z < 4,r H i oc (l+z) 3 ' 37±0 ' 20 . This 
result is about a factor of 1.3 smaller than the commonly 
used t determined by Press et al. (1993) at all z. When 
compared with r at z < 1.5 from HST observations, r at z ~ 
is at least a factor of 4 higher than the one extrapolated 
from at z > 1.5. The different evolution in thi occurs at 
z < 1, a bit lower than that suggested from dn/dz. 

3) For iVm = 10 12 ' 5 ~ 15 cm -2 , the differential column 
density distribution function, /(Ahi), can be best fit by 
/(Nhi) oc with (3 ^ 1.5. The f3 values, however, shows 
the dependence on the A^hi range in the fit, mostly due to 
line blending at Ahi < 10 13 ' 5 cm~ 2 . When combined with 
HST observations, the exponent /3 increases as z decreases 
at < z < 4 for A?hi = 10 13 ~ 17 cm" 2 . 

4) The two-point function of the flux does not repre- 
sent a shape of line profiles as previously suggested. Instead, 
when it is compared at different redshifts, it represents a de- 
gree of line blending in the forest. 

5) The step optical correlation function confirms that 
lines with higher opacity (strong lines) are more clustered 
than lines with lower opacity (weak lines) at the velocity of 
< 200 km s" 1 . 

6) The analyses of the Ly/3 forest at z ~ 2.2 arc in 
good agreement with those of the Lya forest. This result 
shows that previous studies on the Lya forest-only have not 
been significantly biased from line blending and saturation 
at least at z < 2.5. Line blending, however, could be more 
problematic at higher z when it becomes more severe, chang- 
ing some of the results from the analysis of higher- A% i Lya 
forest at N H i > 10 14 ' 5 cm" 2 . 
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APPENDIX A: THE LINE LIST OF Q1101-264 

This is an example of a part of the line lists, which are pub- 
lished electronically at the CDS site to cdsarc.u-strasbg.fr 
(130.79.128.5). 
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Table Al. The line list of Q1101-264 



# 


A 


Identification 


2 


b 


a(b) 


log Nu i 


°\log «Hi) 




( A ) 






km s _1 


km s~ 1 


cm - 2 


cm -2 





3233.622 


H I 


1.65995 


38.76 


13.56 


12.291 


0.124 


1 


3234.807 


H I 


1.66092 


43.36 


9.53 


12.554 


0.073 


2 


3237.851 


Fe II 2382.8 


0.35886 


17.06 


7.92 


12.002 


0.207 


3 


3238.092 


Fe II 2382.8 


0.35896 


5.95 


0.23 


13.443 


0.019 


4 


3238.304 


H I 


1.66380 


9.46 


3.81 


12.115 


0.109 


5 


3238.311 


Fe II 2382.8 


0.35906 


7.14 


4.99 


11.806 


0.186 


6 


3238.499 


Fe II 2382.8 


0.35913 


4.89 


0.50 


12.716 


0.021 


7 


3238.673 


Fe II 2382.8 


0.35921 


5.95 


2.25 


11.979 


0.087 


8 


3238.713 


H I 


1.66414 


15.04 


5.94 


12.100 


0.124 


9 


3240.386 


H I 


1.66551 


30.22 


3.03 


12.558 


0.036 


10 


3245.307 


H I 


1.66956 


26.69 


1.47 


12.734 


0.020 


11 


3249.002 


Fe II 1144.9 


1.83770 


5.63 


0.85 


11.480 


0.036 


12 


3249.567 


Fe II 1144.9 


1.83820 


7.99 


2.73 


11.826 


0.153 


13 


3249.705 


H I 


1.67318 


25.80 


13.77 


12.283 


0.265 


14 


3249.706 


Fe II 1144.9 


1.83832 


5.01 


2.79 


11.660 


0.217 


15 


3249.961 


Fe II 1144.9 


1.83854 


5.86 


0.19 


12.836 


0.011 


16 


3250.156 


Fe II 1144.9 


1.83871 


8.48 


0.62 


12.663 


0.020 


17 


3250.212 


H I 


1.67360 


23.72 


1.35 


13.258 


0.028 


18 


3250.369 


Fe II 1144.9 


1.83890 


5.41 


0.12 


13.161 


0.008 


19 


3250.641 


Fe II 1144.9 


1.83914 


4.84 


0.33 


12.500 


0.016 


20 


3250.758 


Fe II 1144.9 


1.83924 


1.04 


0.99 


12.598 


0.564 


21 


3250.861 


Fe II 1144.9 


1.83933 


6.01 


1.61 


12.056 


0.079 


22 


3251.115 


H I 


1.67434 


13.61 


3.34 


12.217 


0.080 


23 


3252.117 


H I 


1.67516 


22.86 


1.07 


12.932 


0.017 


24 


3252.864 


H I 


1.67578 


22.64 


1.06 


12.916 


0.017 


25 


3256.105 


H I 


1.67844 


48.81 


22.46 


12.651 


0.209 


26 


3256.790 


H I 


1.67901 


24.49 


0.64 


13.910 


0.013 


27 


3257.472 


H I 


1.67957 


26.80 


15.91 


12.204 


0.212 


28 


3261.481 


H I 


1.68287 


14.42 


2.96 


12.935 


0.161 


29 


3261.695 


H I 


1.68304 


33.55 


0.60 


13.999 


0.012 


30 


3262.569 


H I 


1.68376 


33.94 


1.72 


13.458 


0.034 


31 


3262.807 


H I 


1.68396 


15.99 


2.29 


12.909 


0.129 


32 


3265.423 


H I 


1.68611 


51.19 


20.10 


12.551 


0.135 


33 


3266.246 


H I 


1.68679 


19.93 


1.99 


13.370 


0.138 


34 


3266.566 


H I 


1.68705 


25.48 


12.42 


12.869 


0.417 


35 


3271.193 


H I 


1.69086 


32.81 


2.32 


12.693 


0.025 


36 


3273.691 


H I 


1.69291 


33.03 


2.71 


12.765 


0.030 


37 


3275.972 


H I 


1.69479 


28.05 


0.55 


13.536 


0.008 


38 


3278.837 


H I 


1.69714 


68.02 


21.70 


12.852 


0.165 


39 


3279.695 


H I 


1.69785 


30.42 


3.45 


13.047 


0.087 


40 


3280.860 


H I 


1.69881 


119.52 


63.51 


12.912 


0.201 


41 


3283.691 


H I 


1.70114 


52.86 


12.22 


12.709 


0.083 


42 


3284.505 


H I 


1.70181 


23.82 


1.58 


12.985 


0.036 


43 


3288.110 


H I 


1.70477 


59.28 


51.22 


12.563 


0.629 


44 


3288.429 


H I 


1.70503 


14.12 


3.50 


12.702 


0.187 


45 


3288.903 


H I 


1.70542 


40.31 


1.74 


13.971 


0.028 


46 


3289.781 


H I 


1.70615 


20.11 


0.67 


13.497 


0.013 


47 


3290.299 


H I 


1.70657 


16.60 


2.30 


12.623 


0.049 


48 


3292.783 


H I 


1.70862 


31.87 


8.69 


12.314 


0.090 


49 


3295.336 


H I 


1.71072 


31.35 


18.35 


12.089 


0.189 


50 


3296.100 


H I 


1.71134 


20.07 


1.39 


12.857 


0.025 


51 


3296.846 


H I 


1.71196 


10.31 


3.02 


12.078 


0.092 


52 


3298.412 


H I 


1.71325 


48.41 


15.02 


12.435 


0.100 


53 


3299.860 


H I 


1.71444 


23.02 


0.70 


13.511 


0.012 


54 


3301.287 


H I 


1.71561 


17.77 


14.04 


12.464 


0.789 


55 


3301.603 


H I 


1.71587 


18.48 


10.93 


12.975 


0.450 


56 


3302.003 


H I 


1.71620 


30.32 


5.02 


13.362 


0.105 


57 


3304.878 


H I 


1.71856 


26.21 


3.68 


12.421 


0.050 


58 


3306.464 


H I 


1.71987 


23.22 


2.82 


12.618 


0.043 


59 


3309.185 


H I 


1.72211 


35.11 


8.17 


12.313 


0.079 
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